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ON  GOODNESS  OF  FIT  TESTS  FOR  THE  WEI BULL  AND  THE 
EXTREME  VALUE  DISTRIBUTION  WITH  ESTIMATED  PARAMETERS 

by 

Mahesh  Chandra 
Nozer  D.  Singpurwalla 
Michael  A.  Stephens 


1.  Introduction 

The  two  parameter  Weibull  distribution  has  found  many  applications 
in  the  biological,  engineering,  and  the  hydrological  sciences.  For 
instance  it  has  been  used  by  Doll  (1971),  to  describe  the  observed  age 
distribution  of  many  human  cancers.  Its  use  for  describing  failures 
of  electrical  and  mechanical  components  is  well  documented  in  the 
engineering  literature,  and  in  a  comprehensive  study, Benson  (1968) 
discusses  its  use  for  analyzing  flood  data. 

In  this  paper,  we  address  ourselves  to  the  problem  of  testing  the 
null  hypothesis  H^:  that  a  given  random  sample  belongs  to  a  Weibull  or  an 
extreme  value  distribution  with  unknown  parameters.  The  test  statistics  will 
be  EDF  statistics ,  i. e those  based  on  the  empirical  distribution  function, 
and  we  present  tables  of  critical  values  for  testing  H^. 

A  foundation  for  developing  the  tables  of  critical  values  is  the 
recent  theory  by  Durbin  (1973)  on  the  weak  convergence  of  an  "empirical" 
stochastic  process.  This  stochastic  process  is  based  on  the  empirical 
distribution  function  and  estimates  of  the  unknown  parameters.  The 
statistics  that  we  discuss  can  be  represented  as  well-behaved  functionals 


of  this  empirical  process. 
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2.  Preliminaries 

The  two-parameter  Weibull  distribution  is  given  by 


P(T<t)  =  1  -  exp 


t>0 


=  0 


,  otherwise; 


(2.1) 


the  scale  parameter  6  and  the  shape  parameter  ft  are  both  assumed  to  be 
positive. 


If  we  make  the  transformation  X  =  -InT,  where  T  has  the  distribu¬ 
tion  (2.1),  then  X  is  said  to  have  the  extreme  value  distribution. 


P(X£x)  =  F(x)  =  exp 


00  <  x  <  rD, 


(2.2) 


whore  a  =  -ln6  and  b  =  are  the  location  and  the  scale  parameters, 
respectively. 


The  tests  that  we  discuss  in  this  paper  are  for  the  extreme  value 
distribution.  To  make  a  test  of  fit  for  the  Weibull  distribution  we 
shall  take  the  negative  of  the  natural  logarithms  of  the  supposed 
Weibull  data.  Thus,  we  wish  to  consider  the  case  of  testing  whether 
the  distribution  of  a  random  sample  X^,  X^,  ...»  X^,  say  F,  is  an 
extreme  value  distribution  with  unknown  location  and  scale  parameters 
a  and  b  ,  respectively.  Specifically,  we  wish  to  test  the  null 
hypothesis 


Hq:  F(x)  =  G(x) 

for  all  x  and  for  some  (a,b)  ,  where  G(*)  is  the  distribution  F(.) 

given  by  (2.2). 

When  a  and  b  are  specified,  the  Hq  is  said  to  be  "simple," 
and  the  test  reduces  to  testing  the  hypothesis  that  the  independent 

-  2  - 


random  variables 


T-4]  0 


Z  .= 
1 


(y) - -  K-(v ))] 


1  <  i  <  n 


have  a  common  uniform  (0,1)  distribution.  The  Kolmogorov-Smi rnov  test 
is  based  on  the  statistic 


/n  sup  |G  (t)-t| 
n 

0<t<l 


(2.3) 


where 


i 


0  ( t )  =  i 
n  n 


n 

I 

i=l 


.(■(¥!■). 


0  <  t  <  1 ,  and 


(2. A) 


>ihere  1(E)  denotes  the  indicator  of  the  event  E.  Under  the  null 
hypothesis,  the  "empirical"  stochastic  process 


W  (t)  =  /n  (G  (t)-t), 
n  n 


0  <  t  <  i  (2.5) 


satisf ies 


W  (t)  W°  in  «&[0,1]  , 
n 


(2.6) 


d  0 

where-  — >  denotes  convergence  in  distribution,  and  W 

the  Gaussian  process  determined  by  E(W^(t))  =  0,  and 
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li(W°(s)W^(t)  )  =  ciin  (s,t)  -  st,  0  s*  t  1  .  ff  10,1)  denotes  the 
space  of  functions  on  [0,1]  which  are  right  continuous  and  have  left- 
hand  limits. 

When  Hq  is  composite,  an  analogous  test  statistic:  and  a  con¬ 
vergence  theorem  are  obtained;  these  are  discussed  below. 


3.  Asymptotic  Results  When  H^is  Composite 

When  a  and  b  are  not  specified  is  composite,  and  we 

shall  use  (a  ,  S  )  ,  the  maximum  likelihood  estimators  of  (a,b) 
n  n 

Following  Stephens  (1977)  we  call  this  situation  Case  3.  (Cases  1 

and  2  refer  to  the  less  important  cases  in  which  only  a  is  unknown 
or  only  b  is  unknown.) 

Let 


n 


1  <  i  <  n. 


and  define 


H  (t)  =  ~ 
n  n 


E  I(G(Y  .)<t), 
i=l 


t  <  1 


(3.1) 


and 


Y  (t)  =  /n  (H  (t)-t),  0  .  t  <  1  .  (3.2) 

n  n 

Then  from  a  theorem  of  Durbin  (1973)  and  the  appropriate  regularity 
conditions,  the  empirical  process  {y^ ( t) ;  0  <  t  <  l}is  such  that 

Yn  Y°  in  At0-1!  . 

where  Y^  is  a  Gaussian  process  determined  by 
K(Y°(t))=  0 

-  A  - 


,  0  <  t  ^  1 
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and 


K(Y^(s)Y^(t) )  =  min(s,c)  -  st  -  1 . 108(slogs) (tlogt ' 

+  .  257 (slogs) (tlogtlog(-logt)) 

+.  257 (slogs Log (-logs) (tlogt) ) 

-  .60793(slogslog(-logs)tlogtlog(-logt))  ,  0<  s. 

The  above  has  also  been  shown  by  Stephens  (1977).  Tl.e  statist 
interest  in  connection  with  are: 

(i)  the  one-sided  Kolmogorov- Smirnov  statistics 

D+  =  sup  Y  (t)  , 

0<t<l 

D~  =  -inf  Y  (t)  , 

o<t<i  n 

( ii)  the  Kolmogorov- Smirnov  statistic 

D  =  max  (D+,  D  )  , 

(iii)  the  Kuiper  statistic 

V  =  D+  +  6"  , 


0.3) 


ic s  of 


(3. A) 

(3.3) 

(3.6) 

(3.7) 
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(iv)  the  Cramer-Von  Mises  statistics 


2  r1  2 

W  =  /  Y 

Ja  n 


(t)dt  , 


(3.8) 


(v)  the  Watson  statistic 


“2  •  fo  Y"lt,dt  '  IX  V"Ct>dt]* 


(3.9) 


and  the  Anderson-Darling  statistic 


f 1  Yn(t) 

Jn  t(l-t) 


dt 


(3.10) 


As  a  consequence  of  the  continuous  mapping  theorem,  the  limit 

+  —  2  2  2 

laws  of  D  ,  D  ,  D,  V,  W  ,  U  ,  and  A  under  HQ,  are.  given  by  the  laws  o£ 

the  random  variables  sup  Y^(t),  -inf  Y  (t),  max  (D  ,  D  ),  /( Y  (t)  dt, 

O'-tvl  -0-t<l 


0  0 
(*  (t)dt)  - 


Y°(t)dt|  J 


and  lim 


l—  t  Q  O 

f  (Y  (t))  dt  , 

/j  t(l-t) 


respect ively . 


4 .  Sampling  Distributions  of  the  Tost  Statistics 

Stephens  (1977)  has  found  from  theoretical  work  the  cumulants 

2  2  2 

of  the  limiting  distributions  of  W  ,  U  ,  and  A  ,  and  has  us<-d  these 
to  approximate  the  distributions.  The  pth  quantiles  of  these  limiting 
distributions,  together  with  a  modification  for  these  variables  when 
the  sample  size  is  finite  are  given  by  Stephens  (1977,  Table  1). 
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The  sampling  distributions  ol  1)  ,  I)  ,  1)  ,  and  V  have  also  been 
obtained  by  Stephens  (in  some  unpublished  work),  using  Monte  Carlo  methods, 
for  samples  of  size  10,  20,  and  50.  Using  these  finite  sample  results, 

Stephens  uses  an  extrapolation  of  the  quantiles  for  finite  n  ,  to  obtain 
the  corresponding  asymptotic  quantiles.  The  smoothed  Monte  Carlo  points 
are  given  under  Case  3,  in  Table  4.0.  Also  given  in  Table  4.0  are  the  quan¬ 
tiles  of  the  distributions  of  D  ,  1)  ,  D  ,  and  V  when  only  the  scale  parameter 
is  unknown  (Case  1),  and  when  only  the  location  parameter  xs  unknown  (Case  2). 
The  points  for  Case  2  were  obtained  by  Monte  Carlo  methods  similar  to  those 
for  Case  3,  but  for  Case  1,  we  include  some  exact  points  given  by  Durbin  (1975). 

Durbin's  paper  is  concerned  with  testing  for  exponentiality  but,  with 
some  rearrangement,  his  points  apply  for  Case  1.  This  is  because,  in  Case  1, 
b  in  the  extreme  value  distribution  is  known;  then  if  we  make  the  transforma¬ 
tion  y  =  exp(-X/b),  it  is  easily  shown  that  the  distribution  of  y  is  the 
exponential  distribution  F(y)  =  l-exp(-Ay),  y  >  0  ,  where  A  is  exp(a/b). 

Thus  the  problem  reduces  to  testing  that  y  has  the  exponential  distribution 

with  A  unknown  (since  a  is  unknown).  Durbin  has  found  the  exact  points 

for  fnD  ,  /nD  and  /nD  for  this  situation,  and  has  given  extensive  tables 

for  n  ^  100  .  Durbin's  exact  points  have  been  used  where  possible  in 

Table  4.0,  Case  1.  Because  the  transformation  y  =  exp(-X/b)  is  monotonicallv 
decreasing,  the  D+  calculated  directly  from  the  X-values  will  equal  D 
calculated  from  the  y-values,  and  D  (X-values)  will  equal  D+  (y-values) ; 
therefore  the  values  in  Table  1,  Case  1,  for  /nD+  are  Durbin's  values  for 
/nD  and  vice  versa.  Values  of  /nV  in  Table  4.0,  Case  1  are  obtained  from 
Monte  Carlo  methods;  asymptotic  values  for  /nD+,  /nD  and  /nD  are  obtained 
by  extrapolating  Durbin's  exact  points,  bearing  in  mind  that  /nD+  and  /nD 
should  have  the  same  asymptotic  distributions. 


7 


T-410 


The  sampling  distributions  of  1)+,  1)  ,  D  ,  and  V  have  also  been 
obtained  by  Stephens  (in  some  unpublished  work),  using  Monte  Carlo  methods, 
for  samples  of  size  10,  20,  and  50.  Using  these  finite  sample  results, 

Stephens  uses  an  extrapolation  of  the  quantiles  for  finite  n  ,  to  obtain 
the  corresponding  asymptotic  quantiles.  The  smoothed  Monte  Carlo  points 
are  given  under  Case  3,  in  Table  4.0.  Also  given  in  Table  4.0  are  the  quan¬ 
tiles  of  the  distributions  of  D+,  D  ,  D  ,  and  V  when  only  the  scale  parameter 
is  unknown  (Case  1),  and  when  only  the  location  parameter  is  unknown  (Case  2). 
The  points  for  Case  2  were  obtained  by  Monte  Carlo  methods  similar  to  those 
for  Case  3,  but  for  Case  1,  we  include  some  exact  points  given  by  Durbin  (1975). 

Durbin's  paper  is  concerned  with  testing  for  exponentiality  but,  with 
some  rearrangement,  his  points  apply  for  Case  1.  This  is  because,  in  Case  1, 
b  in  the  extreme  value  distribution  is  known;  then  if  we  make  the  transforma¬ 
tion  y  =  exp(-X/b),  it  is  easily  shown  that  the  distribution  of  y  is  the 
exponential  distribution  F(y)  =  l-exp(-Ay),  y  >  0  ,  where  A  is  exp(a/b). 

Thus  the  problem  reduces  to  testing  that  y  has  the  exponential  distribution 

with  A  unknown  (since  a  is  unknown).  Durbin  has  found  the  exact  points 

for  /nD+,  /nD  and  /nD  for  this  situation,  and  has  given  extensive  tables 

for  n  <_  100  .  Durbin's  exact  points  have  been  used  where  possible  in 

Table  4.0,  Case  1.  Because  the  transformation  y  =  exp(-X/b)  is  monotonically 
decreasing,  the  D+  calculated  directly  from  the  X-values  will  equal  D 
calculated  from  the  y-values,  and  D  (X-values)  will  equal  D+  (y-values) ; 
therefore  the  values  in  Table  1,  Case  1,  for  /nD+  are  Dutbin's  values  for 
/nD  and  vice  versa.  Values  of  /nV  in  Table  4.0,  Case  1  are  obtained  from 
Monte  Carlo  methods;  asymptotic  values  for  /nD+,  /nD  and  /nD  are  obtained 
by  extrapolating  Durbin's  exact  points,  bearing  in  mind  chat  /nD+  and  /nD 
should  have  the  same  asymptotic  distributions. 
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A  second  method  of  oht a i n i ng  the  asymptotic  sampling  distributions  of 
tlte  variables  (5.4)  through  (3.10)  is  suggested  in  Wood  (1978),  and  involves 
a  direct  simulation  of  the  Oaussiun  process  Y^.  Specifically,  t  tie  process 
Y^  is  approximated  by  its  finite  dimensional  distribution,  corre¬ 
sponding  to  an  evaluation  of  the  process  at  k  equally  spaced  points  in  the 
unit  interval.  Ten  thousand  multivariate  normal  random  vectors  with 
the  covariance  matrix  given  hy  (3.3)  were  generated  using  the  extended 
precision  version  of  a  program  from  the  IMS  Library.  The  empirical 
distributions  of  the  supremum,  the  inlimum,  and  the  difference  between 
the  supremum  and  the  infimum  of  the  resulting  multivariate  normal  vectors 

wore  then  tabulated,  thus  approximating  the  limit  lavs  of  >n  D+, 

,  _  r.  2  0  1 

>n  D  ,  /n  I),  and  ui  V.  The  Limit  laws  ol  W  ,  U  ,  and  A  were  approximated 
by  using  numerical  integration  techniques.  For  this,  we  used  subroutine 
QSF  from  the  IBM  Scientific  Subroutine  Package.  In  order  to  obtain 
the  quantiles  of  the  true  approximating  limiting  distribution,  (i.e., 
for  for  k  =  00 ) ,  extrapolations  from  finite  values  of  k  must  be 
performed;  how  this  is  done  is  explained  in  Section  4.1. 

Since  Stephens  (1977)  has  already  obtained  the  quantiles  of 

2  2  2. 

the  limiting  distributions  of  W  ,  U'-,  and  A  using  theoretical  methods, 
the  main  purpose  served  by  simulating  the  process  V1  is  to  obtain 
the  quantiles  of  the  limiting  distributions  of  >'n  1)  ,  /n  D  ,  /n  p,  and 
/n  V  by  this  alternative  method.  The  Kolmogorov-Smi rnov  statistics 
D+,  1)  ,  and  1)  are  known  in  similar  goodnoss-ot-f it  situations  to  have 
relatively  low  power.  However,  they  are  commonly  used  in  practice  and 
I)  and  D  are  very  usclul  for  one-sided  tests;  thus  a  comparison  of  the 
two  methods  of  obtaining  quantiles,  both  involving  extropolat ion ,  would 
he  valuable.  We  now  proceed  to  this  comparison. 
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4.1  Results  of  tlie _Monte  Carlo  Simulation  and  extrapolations 

In  ') allies  4.1  througn  4.4,  we  give  the  quant *les  of  the  limiting 
distributions  oi  Lhe  Lest  statistics  (3.4)  througli  (i.10),  obtained 
t rom  statistics  L0,000  replications  of  the  process  Y1'  ,  and  using 
k  =  29,  49,  89  and  119  equally  spaced  points. 

in  order  to  obtain  the  quantiles  of  the  distributions  when  k 

J 

is  infinite,  we  shall  plot  the  kth  quantile  versus  ,  lor  each  of 

V  k 

the  test  statistics  and  extrapolate  to  zero.  The  quantiles  considered 
are  for  p  =  0.75,  0.90,  0.95,  0.975  and  0.99. 

For  example,  in  Figure  4.1  we  show  a  plot  of  the  0.75th  quantile 

versus  — —  ,  for  k  =  29,  59,  89  and  119,  tor  the  test  statistic  l)+  . 
/  k 

The  dotted  line  is  our  linear  extrapolation  to  obtain  the  0.75th 

quantile  of  the  true  limiting  distr ibut ion  (i.e.,  when  k  =  In 

Figure  4.2,  we  show  the  plot  for  the  0.90th  quantile.  Also  shown 

on  the  vertical  axis  of  Figure  4.2  by  an  asterisk,  is  the  0.90th 

quantile  of  tlie  asymptotic  distribution  of  l>+  obtained  by  Stephens, 

and  given  in  Table  4.0.  In  Figures  (4.3)  through  (4.5),  we  show 

analogous  plots  for  p  =  0.95,  0.975,  and  0.99.  Similar  plots,  for 

2  2  2 

the  other  statistics  D  ,  1),  V,  W  ,  11  ,  and  A  ,  ar>'  given  in 
Figures  (4.6)  through  (4.23).  The  asterisks  on  the  vertical  axis  of 
Figures  (4.17)  through  (4.23)  represent  the  appropriate  quantiles  of 
the  limiting  distributions  of  W^,  ligand  A^  obtained  by  Stephens 
(1977,  Table  1).  Since  these  have  been  obLainod  theoretically,  they 
provide  us  with  a  benchmark  for  assessing  the  accuracy  ol  the  simula¬ 
tions  of  tin*  asymptotic  process,  and  also  give  us  some  guidelines  lor 
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Examination  of  Figures  (4.1)  through  (A. 23)  suggests  the 
following  comments: 

(a)  The  plots  of  the  pth  quantile  versus  — —  ,  for 

/k 

all  the  variables  considered  here,  are  approximately 
Linear  for  smaLl  values  of  p  ,  say  p  =  0.75  ,  but 
tend  to  curve  down  for  the  larger  values  ol  p  , 
especially  for  p  =  0.99.  Figure  (A.l)  is  an 
example  of  the  former  and  Figure  4.12  is  an  example 

of  the  latter.  We  should  expect  the  curve  to  be 
.  •  1  ,  , 

monotonic  in  -  ,  and  the  curving  down 

/T 

suggests  that,  as  k  becomes  larger,  the 
accuracy  of  the  simulation  of  Y^  may 
become  suspect . 

(b)  The  linearly  extrapolated  values  shown  in  the 

plots  are,  in  most  instances,  larger  than  the 

corresponding  values  obtained  by  Stephens  (i.e., 

those  indicated  by  the  asterisks).  However, 

parabolic  extrapolation,  also  shown  in  the  plots, 

gives  asymptotic  values  much  closer  to  the 

2  2 

asterisks;  since  the  asterisks  for  W  ,  U  ,  and 
2 

A  can  be  regarded  as  quite  accurate,  it  appears 
that  parabolic  extrapolation  is  to  be  preferred 
to  linear  extrapolation. 
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4 . 2  Simu  La Ling  t  jie_ Maximum  of  a  Brown ian  Motion  Process 


In  comment  (a)  above,  we  have  noted  that  it  appears  that  the 
direct  method  of  simulating  the  asymptotic  process  by  its  discrete 
analogue  evaluated  at.  k  points  may  lead  to  inaccuracies  when  k 
becomes  large.  In  order  to  investigate  this  possibility  further,  it 
was  decided  to  simulate  the  following  Brownian  motion  process,  for 

which  the  distribution  of  the  maximum  is  well  known.  The  process 

n. 

(W(t);  0  <  t  1}  has  mean  0  and  covariance. 

{.E  W(s)W(t) I  =  min  (s,t)  .  (4*1) 


In  order  to  approximate  the  process,  we  obtain  W(t)  at  k  +  1 

equally  spaced  points  on  the  interval  (0,1),  i.e.,  we  evaluate 
% 

lW(j)/k);  j  =0,  1,  ...  k}  by  generating  multivariate  norma]  vectors 
as  described  above.  Siegmund (197  8)  has  shown  Lhat  the  distribution 
of  tiie  maximum  of  this  discretized  process  can  be  approximated  by 


where 


I’{  max  W(j/k)  >  x)  =  2{1  -  d>  (x  + -=4jp-) } 
0  <j  <  k 


0(x) 


1 

/2tT 


Siegmund 's  result 


(4.2) 

is  exact  when  k 


is  infinite. 


In  Table  4.5  we  show  the  quantiles  of  the  distribution  of 
Mk  -  max.  W(/;/k)  ,  0  _  j  v  k  ,  for  k  =  20,  30,  50,  60,  and  90,  obtaiued 
by  simulation,  using  10,000  replications,  and  these  are  compared  with 
tin'  results  given  by  Siegmund’s  approximation  (4.2).  In  Figure  4.24 
these  results  are  shown  graphically.  The  results  of  this  simulation 
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indicate  that  it  would  be  difficult  to  use  the  simulated  points 
aLone  to  obtain  the  correct  intercept  on  the  y-axis;  the  turning- 
down  effect  is  again  present  when  k  becomes  large.  Also,  it 
appears  as  though  either  linear  or  parabolic  ext raoolation  would 
give  reasonable  results,  the  variability  in  the  Monte  Carlo  points 
making  it  difficult  to  distinguish  between  these  methods. 

5.  Quantiles  of  Limiting  Distributions  of  EPF  Statistics 

In  view  of  comment  (b)  above,  supported  by  the  above  results, 
it  seems  reasonable  to  extrapolate  parabolically  to  obtain  asymptotic 
percentage  points  in  Figures  4.1  to  4.16  for  statistics  D+,  D  ,  D  , 
and  V.  The  results  are  given  in  Table  5.1,  for  D+,  D  ,  D  ,  and  V; 
the  values  given  by  Stephens'  Monte  Carlo  method  are  included  for  com¬ 
parison.  It  is  clear  that  there  is  negligible  difference  between  the 
values,  in  terms  of  the  percentage  level,  so  that  the  mean  of  the  two 
estimates  (or  of  all  four  for  D+  and  D  which  should  have  the  same 
quantiles)  might  be  taken  as  a  reasonable  compromise  till  more  accur¬ 
ate  methods  of  finding  true  values  are  available.  Table  5.2  lists 
the  quantiles  of  the  asymptotic  distributions  of  all  the  EDF  statistics 

using  this  compromise  estimate  for  D+,  D  ,  D,  and  V  ,  and  Stephens' 

2  2  2 

theoretically  calculated  values  for  W  ,  U  ,  and  A  . 

6.  Further  Remarks 

(a)  The  above  extensive  study  was  motivated  by  the  desire  to 
compare  two  methods  of  obtaining  asymptotic  quantiles  of  test  statistic 
which  are  functionals  of  a  process  which  is  asymptotically  Gaussian. 
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By  a  comparison  of  the  two  methods,  even  though  each  one  is  based  on 
Monte  Carlo  simulations  followed  by  extrapolation,  points  are  obtained 
which  we  feel  will  be  accurate  for  practical  purposes. 

(b)  The  second  method,  of  directly  approximating  the  asymptotic 
process  by  simulating  a  discretized  version  at  k  points,  is  a  naturally 
appealing  one.  However,  indications  are  that  it  is  very  difficult  to 
preserve  accuracy  as  k  becomes  large.  We  have  to  be  cautious,  since 
we  do  not  always  know  what  to  expect  of  the  calculated  quantity  (in  our 
case,  values  of  functionals  of  the  process)  as  k  becomes  larger.  How¬ 
ever,  we  had  one  good  indication,  given  by  Siegmund's  approximation  for 
the  quantiles  of  the  maximum  of  the  Brownian  motion  process,  which  sug¬ 
gests  that  these  quantiles  should  vary  monotonical ly  in  k  .  This  was 
not  tile  case  for  the  simulated  results;  and  although  we  must  remember 
that  these  are  subject  to  sampling  variations,  the  evidence  overall  in 
these  studies  suggests  that  increasing  k  will  not  necessarily  give 
better  asymptotic  results,  probably  because  the  handling  of  k 
multivariate  normal  vectors  produces  inaccuracies. 
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Table  4.0 

Upper  Tail  Quantiles  for  EDF  Statistics  JxT  Q+,  /a-  D 

/n~  D  ,  i/n  V  .  (For  explanation  of  Cases  1,  2  and  3,  and  the 

Monte  Carlo  technique  Used  see  Section  4) 


Stut  3  Stic 


r~ —  4" 

v  i\  D 

10 

.872 

Case  1 

20 

.878 

50 

.882 

CO 

.886 

/n  D+ 

10 

.988 

Case  2 

20 

1.003 

50 

1.012 

CO 

1.019 

/n  U ' 

10 

.685 

Case  3 

20 

.710 

50 

.111 

(JO 

.12,2 

/n  U 

10 

.773 

Case  1 

20 

.810 

50 

.840 

iXi 

.  8S6 

/ir 

10 

1.012 

Case  2 

20 

1.006 

50 

1.001 

O) 

1.019 

v  n-  0 

10 

.700 

Case  3 

10 

.715 

50 

.724 

oo 

.73 

.969 

1.061 

1.152 

.979 

1.068 

1.176 

.987 

1.070 

1.193 

.994 

1.104 

1.207 

1.135 

1.273 

1.419 

1.152 

1.282 

1.432 

1.168 

1.287 

1.439 

1.174 

1.289 

1.444 

.755 

.842 

.897 

.780 

.859 

.926 

.796 

.870 

.940 

.808 

.876 

.951 

.883 

.987 

1.103 

.921 

1.013 

1.142 

.950 

1.031 

1.171 

.994 

1.104 

1.207 

1.162 

1.275 

1.409 

1,150 

1.280 

1.432 

1.142 

1.290 

1.448 

1.17 

1.296 

1.456 

.766 

.814 

.892 

.785 

.843 

.926 

.796 

.860 

.944 

.81 

.87 

.96 
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Table  4.0  (Continued) 


■  tat  is tie 

u_iP 

.90 

irf  ~  o 

10 

.934 

Case  i 

20 

.954 

50 

.970 

Ou 

.990 

v  n  D 

10 

1.14]. 

Case  2 

20 

1.152 

50 

1.157 

to 

1.161 

/a  i) 

10 

.760 

Case  3 

20 

.779 

50 

.790 

UO 

.797 

v'rT  V 

10 

1.428 

Case  1 

20 

1.460 

50 

1.480 

<Xi 

1.53 

/n  V 

10 

1.388 

Case  2 

20 

1.424 

50 

1.445 

00 

1.459 

/  n  V 

10 

1.287 

Case  3 

20 

1.323 

50 

1.344 

CJ 

1.360 

95 _ .975  .99 


1.026 

1.113 

1.206 

1.049 

1.134 

1.239 

1 . 067 

1.148 

1.263 

1.086 

1.200 

1.300 

1.270 

1.390 

1.520 

1.281 

1.403 

1.525 

1.286 

1.411 

1.528 

1.290 

1.417 

1.530 

.819 

.880 

.944 

.843 

.907 

.973 

.356 

.922 

.988 

.868 

.932 

1.001 

1.547 

1.650 

1.772 

1.575 

1.685 

1.813 

1.593 

1.716 

1.838 

1.65 

1.77 

1.91 

1.493 

1.596 

1.715 

1.538 

1.641 

1.763 

1.564 

1.667 

1.793 

1.584 

1.686 

1.812 

1.381 

1.459 

1.535 

1.428 

1.509 

1.600 

1.453 

1.538 

1.639 

1.471 

1.558 

1.664 
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Table  4.1 


Asymptotic  Quantises*  (Estimates)  of  EDF  Test 
Statxstica  for  the  hx creme  Value  Distribution 
With  Both  Parameters  Estimated 


pth  Quantile 


I  + 

;  j  o 

.  _  .  _L  _ 

D 

D  ; 

V  i 

_ 

wz  1 

t 

_ 

az  ! 

-1 

1 

I  0.010 

1  0.198 

0.197 

0.264  1 

0.489 

0.015 

0.014 

1 

0.098  | 

i 

:  0.025 

i  0.224 

0.226 

0.299 

0.537 

0.018 

0.017 

1 

0.116  | 

'  0.050 

;  0.250 

0.251 

0.326 

0.581 

0.021 

0.020 

1 

0.134  ; 

1 

i  0. 100 

0.280 

i 

0.281 

0.358 

0.637 

0.025 

0.024 

1 

0.158  j 

:  O.230 

!  0.344 

0.343 

0.421 

0.739 

0.035 

0.032 

1 

0.213  | 

i 

i  0.500 

|  0.429 

0.427 

0.503 

0.871 

0.049 

0.047 

0.295 

j  0.730 

1 

'  0.536 

0.532 

0.611 

1.032 

0.072 

0.068 

0.417  ; 

f 

|  0.900 

|  O.ojO 

0.649 

0.  721 

1.199 

0.100 

0.095 

0.569 

1 

j  0.950 

j  0.728 

0.723 

0.791 

1.310 

0.123 

0.117 

0.696 

!  0.975 

j 

j  0.790 

0.794 

0.853 

1.401 

0.146 

0.138 

0.820 

|  0.990 

1 

!  0.868 

0.877 

0.925 

_ 

1.521 

_ 

0.175 

— 

0.165 

0.979 

-’•■Based  on  direct  simulation  of  the  asymptotic  process, 
using  k  =  30  intervals  and  10,000  replicates 
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Table  4.2 

Asymptotic  Quantiles*  (Estimates)  of  EDF  Test 
Statistics  for  the  Extreme  Value  Distribution 
With  Both  Parameters  Estimated 


oth  Quantile 


0.010  j 0.230  i  0.232  0.300  | 0.554 

.256 


j  0.352 

0.636 

0.022 

0.021 

i  0.387 

i 

0.694 

0.026 

0.025 

1 

| 0.452 

0.799 

0.035 

0.034 

!  0.541 

i 

0.941 

0.050 

0.048 

I  0.644 

| 

1 

1.103 

0.073 

0.070 

1  0.753 

1.265 

0.103 

0.097 

I  0.950  |  0.760  j 0.755 

!  0.975  |  0.832  j  0.824 

;  j 

0.990  !  0.916  I  0.923 


0.150 
0.173 

.229  ' 

| 

.316  I 
I 

.446  i 


“Based  on  direct  simulation  of  the  asymptotic  process 
using  k  =  60  intervals  and  10,000  replicates 
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Table  4.3 


Asymptotic  Quantiles*  (Estimates)  of  EDF  Test 
Statistics  for  tae  Extreme  Value  Distribution 
With  Both  Parameters  Estimated 


pth  Quantile 


:  p 

i  =+ 

iT 

_ 

D 

V 

2 

W 

u2 

A2  j 

1  0.010 

' 0.248 

j 

0.244 

0.317 

0.579 

0.016 

0.015 

0.115  j 

!  0.025 

i  0.269 

0.270 

0.342 

0.626 

0.019 

0.018 

0.134 

j  0.050 

!  0.292 

: 

0.295 

0.368 

0.669 

0.022 

0.021 

0.152  : 

j 

i  0.100 

t 

i  0.323 

I 

0.325 

0.400 

0.722 

0.026 

0.025 

0.179  ! 

(  0.250 

j  °-  'ia7 

0.387 

0.463 

0.824 

0.035 

0.034 

0.234  i 

|  0.500 

jO-472 

0.470 

0.551 

0.961 

0.049 

0.047 

0.320 

|  0.750 

j 

;  0.581 

1 

0.579 

0.654 

1.122 

0.073 

0.069 

i 

0,448 

j  0.900 

I  0.694 

0.691 

0.761 

1.293 

0.101 

0.096 

0.610 

'  0.950 

!  0.768 

; 

0.751 

0.828 

1.388 

0.123 

0.117 

0.726 

:  0.975 

j  0.834 

0.823 

0.888 

1.476 

0.143 

0.136 

0.833 

'  0.990 

I 

J 

I - 

j  0.918 

i 

4 - 

0.909 

— 

0.983 

— 

1.597 

— 

0.171 

— 

0.164 

— 

1.010 

'■'Based  on  direct  simulation  of  the  asymptotic  process, 
using  k  =  90  intervals  and  10,000  replicates 
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Table  4.4 


Asymptotic  Quantiles*  (Estimates)  of  EDF  Test 
Statistics  for  the  Extreme  Value  Distribution 
With  Doth  Parameters  Estimated 


! 

| 

pth  Quantile 

1 

1 

| 

P 

i 

*+*  » 

i  !  '  ,2 

•> 

- 1 

2  1 

| 

j 

D  i 

1 

D 

'  D  V  1  W 

1 _ 1 _ _ _ _ 

A 

j 

j  0.010 

i 

!  0.234 

j  0.254 

j  0.328 

J 

;  0.598 

0.017  j 

O.Glo  - 

l 

0.123 

!  0.025 

! 

|  0.2B1 

!  0.278 

j 

!  0.352 

| 

|  0.b38 
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*ilasea  on  direct  simulation  of  the  asymptotic  process, 
using  k  =  120  intervals  an!  10,000  replicates 


-  44  - 


i  vi-  i.  »_  l \.i A  Ci  if.  b 

BENSON ,  M.  A.  (1968).  Uniform  Flood  Frequency  i.stimating  Methods 
for  Federal  Agencies,  Water  Resour.  Res.,  4 (S) ,  891-908. 

JOi.o,  R.  (1971).  The  age  distribution  of  cancer:  implication© 
for  models  of  carcinogens.  J.  Roy.  Statist.  hoc.  her.  A 
133-136. 

DURBIN,  J.  (1973).  Weak  convergence  of  the  sample  distribution 
function  when  parameters  are  estimated.  Ann.  Statist,  f. 
279-290. 

IXkiSIN,  J.  (1975).  Komogorov- Smirnov  tests  when  parameters  are 
estimated  with  applications  to  tests  for  exponentiality  and 
tests  on  spacings.  Biometrika  pj^  5-22. 

SlhCMUND,  D.  (1978).  Corrected  diffusion  approximation  in  certain 
random  walk  problems.  Technical  report  No.  4,  Department  of 
Statistics,  Stanford  University. 

STEPHENS ,  M.  A.  (1977).  Goodness-of-f it  for  the  extreme  value 
distribution.  Biomat rika  pfr  583-588. 

WOOD,  C.  L.  (1978).  A  large-sample  Kolmogorov-Smirnov  test  for 
normality  of  experimental  error  in  randomized  block  design. 
iiiomet  rika  673-676. 


4b 


OFFICE  OF  NAVAL  RESEARCH 
STATISTICS  AND  PROBABILITY  PROGRAM 


BASIC  DISTRIBUTION  LIST 
FOR 

UNCLASSIFIED  TECHNICAL  REPORTS 
JANUARY  1981 


Copies 

Statistics  and  Probability 
Program  (Code  436) 

Office  of  Naval  Research 
Arlington,  VA  22217  3 

Defense  Technical  Information 
Center 

Cameron  Station 

Alexandria,  VA  22314  12 

Office  of  Naval  Research 

New  York  Area  Office 

715  Broadway  -  5th  Floor 

New  York,  New  York  10003  1 

Commanding  Officer 
Office  of  Naval  Research  Eastern/ 
Central  Regional  Office 
Attn:  Director  for  Science 
666  Summer  Street 

Boston,  MA  02210  1 

Commanding  Officer 
Office  of  Naval  Research  Western 
Regional  Office 
Attn:  Dr.  Richard  Lau 
1030  East  Green  Street 
Pasadena,  CA  91101  1 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 

Attn:  Director  for  Science 
536  South  Clark  Street 
Chicago,  Illinois  60605 


Copies 

Office  of  Naval  Research 

Scientific  Liaison  Group 

Attn:  Scientific  Director 

American  Embassy  -  Tokyo 

APO  San  Francisco  96503  1 

Applied  Mathematics  Laboratory 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Attn:  Mr.  G.  H.  Gleissner 
Bethesda,  Maryland  20084  1 

Commandant  of  the  Marine  Corps 
(Code  AX) 

Attn:  Dr.  A.  L.  Slafkosky 
Scientific  Advisor 

Washington,  DC  20380  1 

Director 

National  Security  Agency 
Attn:  Mr.  Stahly  and  Dr.  Maar 
(R31) 

Fort  Meade,  MD  20755  2 

Navy  Library 

National  Space  Technology  Laboratory 

Attn:  Navy  Librarian 

Bay  St.  Louis,  MS  39522  1 

U.  S.  Army  Research  Office 
P.0.  Box  12211 
Attn:  Dr.  J.  Chandra 
Research  Triangle  Park,  NC 
27706  1 


1 


Copies 


Copies 


OASD  ( I&L) ,  Pentagon 

Attn:  Mr.  Charles  S.  Smith 

Washington,  DC  20301  1 

ARI  Field  Unit-USAREUR 
Attn:  Library 

c/o  ODCSPER 

HQ  USAEREUR  &  7th  Army 

APO  New  York  09403  1 

Naval  Underwater  Systems  Center 
Attn:  Dr.  Derrill  J.  Bordelon 
Code  21 

Newport,  Rhode  Island  02840  1 

Library,  Code  1424 

Naval  Postgraduate  School 

Monterey,  CA  93940  1 

Technical  Information  Division 
Naval  Research  Laboratory 
Washington,  DC  20375  1 

Dr.  Barbara  Bailar 
Associate  Director,  Statistical 
Standards 
Bureau  of  Census 

Washington,  DC  20233  1 

Director 

AMSAA 

Attn:  DRXSY-MP,  H.  Cohen 
Aberdeen  Proving  Ground,  MD 
21005  1 

Dr .  Gerhard  Heiche 
Naval  Air  Systems  Command 
(NAIR  03) 

Jefferson  Plaza  No.  1 
Arlington,  VA  20360  1 

B.  E.  Clark 
RR  it 2,  Box  647-B 

Graham,  NC  27253  1 

Leon  Slavin 

Naval  Sea  Systems  Command 
(NSEA  05H) 

Crystal  Mall  #4,  Rm.  129 
Washington,  DC  20036  1 


ATAA-SL,  Library 

U.  S.  Army  TRAD0C  Systems  Analysis 
Activity 

Department  of  the  Army 
White  Sands  Missile  Range,  NM 
88002  1 

Technical  Library 

Naval  Ordnance  Station 

Indian  Head,  MD  20640  1 

Defense  Logistics  Studies 
Information  Exchange 
Army  Logistics  Management  Center 
Attn:  Mr.  J.  Dowling 
Fort  Lee,  VA  23801  1 

Reliabil ity  Analysis  Center  (RAC) 

RADC/RBRAC 

Attn:  I.L.  Krulac 

Data  Coordinator 
Government  Programs 

Griffiss  AFB,  New  York  13441  1 

Mr.  F.  R.  Del  Priori 
Code  224 

Operational  Test  and  Evaluation 
Force  (OPTEVFOR) 

Norfolk,  VA  23511  1 

Mr.  Jim  Cates 
Code  9211 

Fleet  Material  Support  Office 
U.  S.  Navy  Supply  Center 
Mechanisburg,  PA  17055  1 

Mr.  Ted  Tupper 
Code  M-311C 

Military  Sealift  Command 

Department  of  the  Navy 

Washington,  DC  20390  1 

Mr.  Barnard  H.  Bissinger 
Mathematical  Sciences 
Capitol  Campus 

Pennsylvania  State  University 
Middletown,  PA  17057  1 


Morton  L.  Buckberg 
Naval  Sea  Systems  Command 
(NSEA  3133) 

National  Center  No.  2 

Arlington,  VA  ?0360  1 


-  2  - 


